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ABSTRACT 



We use the 2 Micron All Sky Survey (2MASS) extended source catalogue to de- 
termine area-averaged angular correlation functions, Wp, to high orders {p ^9). The 
main sample used contains 650 745 galaxies below an extinction-corrected magnitude 
of Ks=H.b and limited to |6| >10° and represents an order of magnitude increase in 
solid angle over previous samples used in such analysis (Ri34 000 deg^). The high-order 
correlation functions are used to determine the projected and real space hierarchi- 
cal amplitudes, Sp = u}p/u}2~^ and Sp = Cp/'ff"^- I'^ contrast to recent results, for 
j3 ^ 6 these parameters are found to be quite constant over a wide range of scales 
to r «40/i~^Mpc, consistent with a Gaussian form to the primordial distribution of 
density fluctuations which has evolved under the action of gravitational instability. We 
test the sensitivity of our results to the presence of rare fluctuations in the local galaxy 
distribution by cutting various regions of over-density from the main sample; unlike 
previous analyses, we find that at least for p< 4, the high-order clustering statistics are 
relatively robust to the removal of the largest superclusters. Since we probe well into 
the linear regime we are able to make direct comparisons with perturbation theory; we 
use our constraints on the Kg-hand Sp parameters in two ways. First, we examine their 
consistency with non-Gaussian initial conditions; we are able to rule out strong non- 
Gaussianity in the primordial density field, as might be seeded by topological defects 

— — r)/2 

such as cosmic strings or global textures, producing relations of the form ^p oc ^2 i 
the ?s2.5ct confidence level. Second, we investigate the way in which galaxies trace the 
underlying mass distribution. We find evidence for a non-zero quadratic contribution 
to the galaxy bias, parameterised by C2=0.57±0.33. This positive result represents a 
significant difference from the negative values found previously; we examine a possi- 
ble explanation in the light of recent observations which universally provide negative 
values for C2. 

Key words: surveys - galaxies: clustering - cosmology: observations - large-scale 
structure of the Universe - infrared: galaxies 



1 INTRODUCTION 

The variance of local galaxy density fluctuations has been 
measured to high accuracy over the last decade, both 
throu g h the 2-point corre la tion function ( e .g. [ Zehavi et alJ 
I2OO4I: iMaUer et aP l2005l : iHawkins et"ai] l200 3h and its 
Fourier transform, the power spectrum fe.g. lCole ct al. 2005; 
iFrith et alJl20 05m. For a perfectly Gaussian density field, 
the 2-point statistic forms a complete description of the 
galaxy distribution as all higher order connected moments 
are zero. 
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Assuming a Gaussian form for the primordial distri- 
bution of density fluctuations, perturbation theory predicts 
non-zero higher order correlation functions of the local 
galaxy density fleld; as the initial inhomogeneities grow grav- 
itationally, the density distribution becomes asymmetric, de- 
veloping non-zero skewness and kurtosis etc. However, it is 
possible to construct models of inflation, and also models 
which contain non-linear structures in the primordial den- 
sity field, such that the initial conditions themselves are non- 
Gaussian (e.g. _Silk & Juszkiowicz 1991; Weinberg & Col3 
Il992ll . Therefore, if non-zero high-order moments of the lo- 
cal galaxy density field are detected, these arise either via 
the gravitational collapse of initially Gaussian density fluc- 
tuations or from some degree of primordial non-Gaussianity. 



2 W.J. Frith, P.J. Outram & T. Shanks 



Examining the high-order correlation functions of the 
local galaxy distribution therefore represents a particularly 
powerful approach to understanding the nature of the pri- 
mordial density field. It is now well established that under 
the assumption of Gaussian initial conditions, these high- 
order moments are expected to display a hierarchical scal- 
ing such that the p-order cumulants of the local density 
field < 5^ >c can be expressed in terms of the variance of 
the distributi o n, such that < 5^ >e= < S ' ^ >''~^ (e.g 
i Peebled llQSd: IfVvI|1984I: iBouchet et al.J fl992l: iBernardeaul 
Il99lll994al: l( 
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These Sp coefficients, which quantify the departure from 
Gaussian behaviour, are insensitive to cosmic epoch or as- 
sumed cosmology at scales where the growth of the den- 
sity field is linear or quasi-linear. Departures from the hi- 
erarchical scaling of the high-order galaxy correlation func- 
tions on these scales might be expected only in the case 
of strongly non-Gaussian initial conditions or some form 
of scale-dependency at large scales in the bias between the 
galaxy distribution and the underlying mass distribution. 

Several surveys have made a significant contribution to 
our understanding of this issue. It has long been known 
that high-or der correlation function s are non-zero on small 
scales (e.g. iGroth fc Peebles! Il977l: I Saunde rs et al.l [l991 : 



Bouchet et alj Il993l: iGaztaiiagal Il994: Szapu di et all boOC ; 
Hovle et alJl2000t) . More recentlv Tcroton et al. " ( 2004") used 



used 

volume-limited catalogues drawn from the 2dF Galaxy Red- 
shift Survey (2dFGRS) to constrain Sp in redshift space 
(for p ^5) to r RilO/i"^Mpc. For r<4/i"^Mpc the hier- 
archical scaling parameters, Sp, were found to be approxi- 
mately constant. However, at larger values of r these coeffi- 
cients increased with scale for some luminosities, consistent 
with some models of structure formation with strongly non- 
Gaussian initial density fields ("Gaztanaga & Mahonen^^9d 
[Gazta fiaga & Fosalba 1998; White 1999; Bcrnard eau et alJ 
120021) . However, the results were complicated by the pres- 
ence of two massive superclusters which, when removed from 
the analysis (corresponding to a reduction in the volume of 
~2 per cent), now resulted in constant scaling parameters 
on all scales to r ~W Mpc, consistent with Gaussian ini- 
tial conditions. It appears therefore that the 2dFGRS does 
not probe a large enough volume to constitute a fair sam- 
ple of the local galaxy distribution for high-order correlation 
functions. 

The form of the scaling parame ters on larg e r scal es are 
also poorly constrained. Previously. iGaztanagal il994l) used 
the APM galaxy survey to bj = 20 to constrain high-order 
(p <9) angular correlation functions and the associated an- 
gular scaling parameters, Sp, to d «7° (corresponding to 
r< 30 Mpc). Despite the fact that the APM galaxy sur- 
vey covers a ~4x larger solid angle than the 2dFGRS, a 
similar increase in S3 on large scales was found as seen for 
S3 observed in the 2dFGRS. At higher orders, further de- 
partures from the expected hierarchical scaling were also ob- 
served, although the scales probed were limited to 6< 3° and 
the statistics became increasingly uncertain. These features 
have also been measured in the smaller Edinburgh-Durham 
Galaxy Catalogue (Szapudi fc Gaztanaga 1998) . Of course, 
the in crease in S3 at large scales detected by iGaztafiagal 
(|l994|) may also be due to the same massive supercluster 
observed in the Southern 2dFGRS field, as the 2dFGRS 
is drawn from the larger APM galaxy sample. However, it 



would be surprising if the volume probed by the APM galaxy 
survey still did not constitute a fair sample of the Universe. 

High-order correlation functions and the amplitude of 
the associated scaling parameters also provide us with a 
powerful probe of the way in which galeixies trace the under- 
lying mass distribution. Recent measurements of the linear 
bias, that is the bias associated with the variance of the 
density field such that ^2,gai = fe?^2.DM, indicate that in 
the near infra-red &i >1 (e.g. iFrith et alJll2005b'l whereas 
for optically-selected g alaxies bi (e.g. IVerde et alJl2002l : 
iGaztanaga et al]l2005h . Whether there exist non-linear con- 
tributions to the galaxy bias, such that the bias is a func- 
tion of the density field, can be examine d using high-order 
moments of the galaxy density field (e.g. iFrv fc Gaztanagal 
Il993h . 

Previous high-order clustering analyses have universally 
found that the observed skewness etc. are significantly lower 
than the expected values in a ACDM cosmology. This dis- 
crepancy has been interpreted as evidence for a negative, 
non-linear contribution to the galax;y bias. Most recently, 
'Gaztafiag a et all ll2005l) used volume-limited samples of the 
2dFGRS to determine the redshift space 3-point correlation 
function and constrained bj-band bias parameters to bi = 
0.94^Q jj and a non-linear, quadratic bias of C2 = b2/bi = 
— 0.36lgjg, although as with previous 2dFGRS analyses 
these results are significantly affected by the presence of 
t wo m assive superclusters. Independently, |Pan fc Szapudi 
(||200^ determined the monopole contribution to the 2dF- 
GRS 3-point correlation function and determined a similar 
&,j-band linear bias to Gaztanaga et al. (2005) and a neg- 
ative C2 parameter as well, although considerably smaller 
and at a reduced significance to the other 2dFGRS con- 
straint. Computing the bispectrum for the PSCz catalogue 
(which is selected from the 60/^m IRAS galaxy sample), 
iFeldmaiTet al. ( 200lh also constrain infra-red bias param- 
eters of bi = 0.83 ± 0.13 and C2 = -0.50 ± 0.48. 

The 2 Micron All Sky Survey (2MASS) has recently 
been completed and provides Ks, H and J-band photom- 
etry for 1.6x10® extended sources over the entire sky to 
-R's> 13.5. 2MASS is the largest existing all-sky galaxy sur- 
vey and therefore represents a uniquely powerful probe of 
the local galaxy density field at large scales; the solid angle 
of the 2MASS |6| ^10° sample used in this paper represents 
an order of magnitude increase over the APM galaxy sur- 
vey, meaning that clustering statistics determined from the 
2MASS galaxy sample will suffer less from projection effects 
while probing a comparable volume. The 2MASS survey also 
represents an order of magnitude increase in volume over the 
largest volume-limited 2dFGRS sample; it is therefore pos- 
sible to probe much larger scales than the 2dFGRS although 
with the added complication of projection effects. A further 
advantage of 2MASS over previous datasets is that the pho- 
tometry is extremely accurate with high completeness for 
|6| ^10°; correcting for the variable completeness over the 
2dFGRS survey area for instance, complicates the analysis 
and increases the uncertainty from possible systematic ef- 
fects. The main drawback to high-order clustering analysis 
of the 2MASS data (as with the APM galaxy survey) is the 
lack of availiable three-dimensional information; the clus- 
tering signal from a particular scale in real space is smeared 
over a range of angular scales. For this reason detailed fea- 
tures in the real space correlation function, such as the shoul- 
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der at r »10/t~^ Mpc iBauehl [l993 : iHawkins et alJ 120031: 
IZehavi et alJ2004) . may not be detected by the 2MASS pro- 
jected correlation function. 

In this work, we aim to determine the high-order angu- 
lar correlation functions and the associated scaling param- 
eters of the local galaxy density field to high precision and 
large scales, using the final 2MASS extended source cata- 
logue. In section 2, the details of the method for estimating 
high-order correlation functions are given. The galaxy sam- 
ple used and the error analysis are described in section 3. 
We present the p-point angular correlation functions (for 
p ^9) and the hierarchical scaling relations in section 4. 
We also examine possible systematic eff'ects arising from ex- 
treme fluctuations in the observed galaxy density field. The 
implications for the form of primordial density fluctuations 
and non-linear galaxy bias are discussed in section 5. The 
conclusions follow in section 6. 



2 METHOD OF ESTIMATION 

2.1 The p-point correlation function 

The p-point galaxy correlation function estimates the joint 
probability that p galaxies are separated by a certain scale 
and can be defined through considering fluctuations in the 
galaxy density field. The connected or reduced part of this 
statistic corresponds to the contribution to this probability 
which does not include any conditional probability on lower 
orders: 



LJp{dl 



, 5„ >c 



(1) 



where 5 denotes the density fluctuation; for p ^3, the unre- 
duced and reduced correlation functions are the same. In 
this paper, we work with the reduced p-point correlation 
function only. 

The 2-point angular galaxy correlation function, CJ2(S), 
is given in terms of the probability of finding two galaxies 
in area elements dfii and dQ2 separated by angle 6: 

dP2 = A/"^ [H- 0^2(61)] dflidQz, (2) 

whe re Af is the mean number of galaxies p er unit solid angle 
fe.g. lGroth fc Peebleslll977l:lPeebleslll980h . Similarly, the 3- 
point function, u!3{9), is defined by the joint probability of 
finding galaxies in each of three area elements: 

dPs = A/'Ml + cJ2iei2) + 0^2(6*23) + C^2(6ll3) 

+ uj3{ei2,023,e3i)] dnidn2dn3 (3) 

The first terms in equations |5| and |3 correspond to the con- 
tributions from galaxy pairs or triplets respectively which 
are accidentally seen as close together in projection but are 
at very different radial distances. Similarly, the following 
three terms in equation |3 account for the contribution from 
one correlated pair and a third uncorrelated galaxy which 
forms a triplet by chance line-of-sight clustering. This leaves 
the final term which defines the contribution from the real 
clustering of triplets. 

A simple way in which to estimate the high-order cor- 
relation functions is through the area-averaged correlation 
function, cop, defined as: 



ijp(e) = — 



LUp{6i,..., 6p) dVli, 



, dQ,^ 



(4) 



where Q. is the solid angle of the cone defined by its angular 
radius 6. The area-averaged correlation function, ujp, is es- 
timated by considering the central moments of the angular 
counts: 



mp{e) =< (iV - iV)" 



> 



N=ao 



(5) 



where Pn(Q) denotes the count probability distribution 
function and is calculated by placing circular cells of angular 
radius 6 over the survey area and determing the number of 
cells containing exactly N galaxies: 

Nn 
Nc' 



PN{e) 



(6) 



where A'jv and Nc denote the number of cells containing 
N galaxies and the total number of cells respectively. N in 
equation|K|is the mean number of galaxies in a cell of angular 
radius 6 and may be determined directly from the count 
probability distribution function: 



N = 



N—oc 

E 

JV=0 



NPN{e) 



(7) 



The moments of the count probability distribution func- 
tion determined via equation |K| yield the unreduced corre- 
lation function through the relation rrip =< 5'' > N'' . In 
order to obtain the reduced correlation function, the con- 
nected moments Hp are determined: 

fJ-2 = m2 
fi3 = ms 



^4 = m4 — 3m2 



(8) 

(see iGaztafiagal il994) for higher order relations). In ad- 
dition we apply a shot noise correction iGaztafiagal 1199^ : 
iBaueh et aI^ll995^ such that: 



k2 


= M2 


-iV 




k3 


= M3 


-3fc2 


-N 


ki 


= M4 


-7k2 


-6k3-N 



(9) 

The area-averaged, reduced angular correlation function is 
then determined from the relation cop —< 5^ >c= kp/N^ . 



2.2 Hierarchical scaling 

In perturbation theory, the density field, evolved by gravity 
from an initially Gaussian distribution, leads to a hierarchi- 
cal clustering pattern, such that all high-order correlations 
can be expressed in terms of the 2-point correlation function: 



(10) 



where ^p,dm is the volume-averaged p-point dark matter 
correlation function. Importantly, this relation is preserved 
for the galaxy density fi eld such that (.p^gai ~ Sp^gai ^290! 
iFrv fc GaztafjagtM Il993l) . The relation between the dark 
matter and galaxy density fields may be expressed through 
a Taylor expansion of the dark matter density contrast: 



5gal = ^ ^(cSdm)'' 



(11) 



p=0 
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For t he skewness, it can be shown that iFrv fc Gaztanag3 
Il993l): 



•S'a.gal 



- (S3, DM + 3C2), 



(12) 



where 61 is the hnear bias such that £,2,gai = -Dm on 
scales where the variance ol the density field is small, and 
C2 = 62/61 quantifies the second-order contribution to the 
galaxy bias. 

Using perturbation theory, it is possible to determine 
precise quantitative predictions for the Sp^oM parameters. 
Assuming a power law form for the three-dimensional power 
spectrum of density fluctuations, P{k) oc fc", and a top-hat 
window function, the skew ness of the matter distr ibution, 
Ss.DM may be determined djuszkiewicz et al]ll993^ : 



c _ 34 

03, DM — — 



(n+3) 



(13) 



We later use these expressions to determine constraints on 
non-hnear bias through comparisons with Ss^gai in the lin- 
ear and quasi-li near regimes. We as s ume a power spe ctrum 
slope of n=-2 jPercival et alJl200ll : |CoieetaO|200i). The 
uncertainties on n are small (An <0.1) compared to the 
sampling errors in the measurement of S^^gai', we therefore 
assume the concordance value cited above to generate pre- 
dictions for S3, DM and neglect small uncertainties in this 
parameter. 



2.3 Transformation to three dimensions 

For a projected galaxy distribution, a hierarchical scaling 
relation may also be defined in terms of the area-averaged 
p-point galaxy correlation functions and angular scaling co- 
efficients, s„: 



-,p-i 



(14) 



We wish to transform these angular scaling parameters, 
Sp, to the three-dimensional coefficients, Sp, in order to 
make comparisons with perturbation theory and constrain 
non- li near galaxy bias. Following the method of Gaztaf iagai 
il994ri. we transform to three dimensions via the relation: 



Spi0)Bp{j) 
'■p(7)Cp(7)' 



(15) 



where f — OT) is the mean scale probed at an angular 
scale 9 for a survey of median depth 2?, 7 is the slope of 
the 2-point real space correlation function, and Bp and Cp 
are related to the number of different configurations of the 
three-dime nsional and angu lar hierarchical tree graphs re- 
spectively (iGaztanagal I1994I) . Here we use the form of the 
real s pace correlation func tion determined from the 2dF- 
GRS (|Hawkins et al.ll2003^ . The rp factor is related to the 
selection function ^Z: 

_ ij^ Ip 



TP~1 



(16) 



where 



>F[x)x^dx (17) 



where e describes the e volution of clustering with redshift 
jGroth fc Peebles! Il977ll and is taken to be e=0 in good 



agreement with recent observational and theoretical con- 
siderations jHamihon et alJll99ll: iPeacock fc Doddj Il994l: 



ICarlberg et alj|2000l : IWilsonI 



2003). Here, x denotes the co- 



moving distance and F{x) a correction for curvature such 
that ^'(a;) = [1 - {Hox/cf{Q,m - 1)]^^^. In this work we use 
the concordance value of Jim =0.3. For the selection function 
we use a parameterised form for the n{z) such that: 



(-(^)"1 



^ 2(.71.412)3 

ijBaugh fc EfstathioxJ Il993l: iMaUer et aP l200,4 l where z is 
determine d from the 2MASS-2 dFGRS matched sample de- 
scribed in iFrith et"al For reference 2=0.074 for 
Ks <13.5. In this case, the normalisation of the selection 
function factors out (see equations 1161 and I17I I. 

The transformation described in equation 1151 is robust 
to reasonable changes in t he selection function and choice of 
cosmological parameters iGaztafiagal [l994'l . However, this 
relation becomes uncertain on large angular scales, 6 >2° 
ySaztanaga 1 993) . due to the fact that the power law form 
to the 2-point correlation function and the value of 7 are 
not well constrained on large scales. 



3 ANALYSIS OF THE 2MASS DATA 

3.1 The 2MASS Extended Source Catalogue 

We select objects from the 2MASS final release extended 
source catalogue (Jarrctt 200^1) above a galactic latitude of 
|6| = 10° in order to remove regions of high extinction and 
stellar contamination (see Fig. |SJ and below an extinction- 
co rrected magnit u de lim it of its = 13.5 using the dust maps 
of ISchlegel et al.l Jl998l) . The 2MASS dataset removes or 
fiags sources identifi ed as artefacts such as diffraction spikes 
and meteor streaks Jjarrett et alJl2000t) : we use the 2MASS 
cc.flag to remove such objects. We also employ a colour 
cut (J - Ks <0.7 and J - Ks >1.4) below Ks = 12 in or- 
der to remove a small number of objects identified as non- 
extragalactic extended sources tMaller et al . 2003, 200^. 
The subsequent sample of 650 745 galaxies probes to a me- 
dian depth of "D ~220 Mpc and covers 83 per cent of the 
entire sky («34 000 deg^). 

We estimate the total iC'a-band fiu x using the magni- 
tude estimato r described inlFrith et al.l ll2005a b). based on 
the method o f ICole et alJj200J) . To recap, we use the deeper 
and more accurate 2MASS J-band extrapolated magnitudes 
and colour-correct to the iC's-band using the J and Ks- 
band fiducial elliptical Kron magnitudes; this provides ex- 
ce llent agreement with the independent K-hand photometry 
of lLovedavl ll2000tl and also the 2M ASS-selected 6dF Galaxy 
Survey total magnitude estimator j Jones et alj|2003) . 

In order to calculate the high-order correlation func- 
tions of this 2MASS sample, we determine the count prob- 
ability distribution function detailed in equation |S| by ran- 
domly placing Nc = W^ cell centres over the survey area. 
Each cell is then allowed to grow and the number of galax- 
ies as a function of the angular radius is recorded. The size of 
the cells is limited by the galactic latitude limit of \b\ = 10° 
in the 2MASS sample. We replace cells which are lost as 
they encroach the boundary of the galaxy sample such that 
A'^c=10^, independent of the angular scale probed. We use 
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cell radii in the range 0°.01 to 25°. 1, equivalent to a range 
in the mean scale probed of 0.04 Mpc to 104 Mpc. 



3.2 Error estimation 

The statistical uncertainty associated with the correlation 
function and angular scaling parameters in this work are de- 
termined using bootstrap estimates. The full 2MASS sample 
is split into 20 equal area regions of «1700 deg^; 20 of these 
sub-areas are selected at random (with repeats) and the as- 
sociated clustering statistic is determined. We repeat this 
1000 times and determine the standard deviation; these are 
indicated by the errorbars in Figs.0and|5| The size of these 
errors is not significantly altered if we vary the number of 
realisations or sub-areas within reasonable limits. 

Once we have determined the angular correlation func- 
tions and scaling parameters and the associated errors in 
this way, we wish to carry out comparisons with predic- 
tions from perturbation theory. Since correlation function 
estimates determined at different cell radii are highly cor- 
related being integral quantities, it is necessary to account 
for the covariance between each datapoint when performing 
fits to the data. We determine the covariance matrix (Cij) 
using the bootstrap method above. We then determine x^- 

x' = ^5]A,CriA,, (19) 

i 3 

where Ai = s'^"{i) — «""''(*) for example. We compute 
the inverse covariance matrix, C~j^ , using the Numerical 
Re cipes Singular Value Decomposition algorithm. As noted 
by ICroton et al.l ll2004fl , previous constraints on the high- 
order scaling parameters (except for the 2dFGRS results) 
ignored the correlations between different bins leading to 
unrealistically small errors in the fitted values. 



p 


Sp 


Sp 


3 


5.28 ± 0.45 


4.29 ± 0.39 


4 


57.7 ± 9.2 


36.0 ± 5.7 


5 


1510 ± 507 


678 ± 228 


6 


(3.68 ± 2.08) X 10"' 


(1.20 ± 0.68)xl0"' 


7 


(9.74 ± 8.72) xIQS 


(2.25 ± 2.01)xl05 


8 


(2.64 ± 3.69) xlO'^ 


(4.19 ± 5.86)xl0'' 


9 


(2.03 ± 1.56) xlO^ 


(2.26 ± 1.74) X 10* 



Table 1. Small scale fits to the high-order scaling parameters 
(assuming constant values) determined for the full 2MASS |6| ^ 
10°, Ks <13.5 sample (see Fig. 0. In the second column the 
best fit angular scaling parameters (for 3^ p ^9) are shown, 
determined from fits in the range 0°.04 < d <1°.0 for p ^8 
and 0°.04 < <0°.25 for p=9. Similarly, in the third column we 
show the best fit real space scaling parameters, Sp, fitted over 
an equivalent range of scales, 0.15< f <4.0/i~^Mpc for p ^8 
and 0.15< f <1.0h~^ Mpc for p=9. The errors in each case take 
into account the covariance matrix determined from the bootstrap 
estimates described in section 3.2. 
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Sp 


Sp 


3 


4.91 ± 0.60 


4.00 ± 0.49 


4 


54.2 ± 11.1 


33.8 ± 6.9 


5 


740 ± 240 


332 ± 108 


6 


(1.01 ± 0.67) X 10"' 


(3.30 ± 2.19)xl03 



Table 2. Large scale fits to the high-order scaling parameters (as- 
suming constant values) determined for the full 2MASS |6| 10°, 
Kg <13.5 sample (see Fig.0. In the second column the best fit 
angular scaling parameters (for 3^ p ^6) are shown, determined 
from fits in the range 1°.0 < 9 <10°. Similarly, in the third 
column we show the best fit real space scaling parameters, Sp, 
fitted over an equivalent range of scales, 4.0< f <40/i~^Mpc. 
The errors are determined as in Table ITl 



4 RESULTS 

4.1 Area-averaged correlation functions 

In the first column of Fig. we plot the area-averaged cor- 
relation functions for p SC9 determined for the full |b| ^10° 
2MASS galaxy sample described in section 3.1. In each case 
the errors are determined via the bootstrap technique de- 
scribed in section 3.2. For reference we compare our result 
for ijj2 with the best f it to the 2MASS K s <13.5 2-point 
correlation function of lMaller et~al] ^2005^ (A=0.10±0.01, 
l-7=-0.79±0.02, determined on scales of 6 ^2°. 5 where 
U2 = AO^-''). 

The p-point correlation functions are consistent with 
power laws over several orders of magnitude in angular 
scale. In particular, we note that the area-averaged ijJ2 deter- 
mined in this work is in good agreemen t with the directly- 
determined result of lMaller et alj ^2005^ at small scales (we 
find A = 0.13 ± 0.01, l-7=-0.77±0.04 for 9 s:2°.5). This 
agreement is robust to ch anges in th e magnitude estimator 
and galactic latitude cut; iMaller et a l. (2005) use the Ks- 
band elliptical Kron magnitude estimator and a |6| ^ 20° 
cut. 

In order to determine whether these results are consis- 
tent with the hierarchical scaling described in section 2.2 



we compute the angular scaling parameters, Sp (see equa- 
tion ^J; these are shown in the second column of Fig. 
The angular scaling parameters are transformed into the 
real space Sp parameters in the third column of Fig. (see 
eguation lldH . The Sp and Sp coefficients are consistent with 
constant values to large scales (^< 20° , r< 100 Mpc), al- 
though there may be a decrease in amplitude between the 
non-linear and quasi-linear regime (6 «1°, r «4/i-iMpc). 
One caveat to this is that while these coefficients are consis- 
tent with a slope of zero, the constraints become increasingly 
weak at higher orders such that in the range 1°.0 < ^ <10° 
the slopes are constrained to 7(53) = O.OI^q ^si li^i) = 
QmtVm, 7(S5) = -O.Oll°;^^ ^(^^^ ^ -0.39t°J« at la 
confidence (where Sp (x 6~' , marginalising over the normali- 
sation) ; we investigate the constraint on the slope of 5*3 and 
the level to which this can reject primordial non-Gaussianity 
in section 5.1. 

Since we probe well into the linear and quasi-linear 
regimes we are able to make comparisons with predictions 
from perturbation theory. We perform fits to these scaling 
parameters considering the covariance in the datapoints (see 
equation 1191 . We consider small and large scales separately 
due to the fact that, despite the consistency of the scaling 
parameters over all scales considered, pert urbation theory is 
not expected to be valid on small scales jBernardeau et all 
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Figure 1. Moments of the 2MASS galaxy density field for the full |6| ^ 10° Ks <13.5 sample (all on logarithmic scales). Each row 
corresponds to a diff'erent moment of the galaxy density field. Since in each column we plot a different statistic {tOp, Sp or Sp) we denote 
the X-axis label as Fp. In each case the errors are taken from bootstrap estimates described in section 3.2. Datapoints with extremely 
large error bars are omitted for clarity. The columns are set out as follows: 

First colu mn: Area-averaged correlation functions for 650 745 Ks <13.5 2MASS galaxies. In the 012 panel, the best fit result at small 
scales from lMaller et ahl l2005^ is indicated by the dotted line. 

Second column: The angular scaling parameters (sp) determined via equation 1141 for the full sample. 
Third column: The real space scaling parameters (Sp) determined via eauation ll5l for the full sample. 
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Figure 2. Moments of the 2MASS galaxy density field making various cuts to the full |fe| ^ 10° Kb <13.5 sample. As in Fig.lTlcach row 
corresponds to a different moment and the errors are determined via bootstrap estimates as described in section 3.2. For p=2 (the top 
row) we show the area-averaged correlation function; for higher order moments we display the angular scaling parameters, Sp. In each 
case we indicate the corresponding results from the full sample shown in Fig. 0by a solid line. Datapoints with extremely large errorbars 
are omitted for clarity. The columns are set out as follows; 

First column: We show the 2-point function and the higher-order angular scaling parameters having removed a region of radius 6° 
centred on the largest over-density in the sample, the Shapley supercluster (see section 4.2). This corresponds to a removal of 1.1 per cent 
of the galaxies and 0.33 per cent of the solid angle of the full sample. 

Second column: We apply a more stringent cut such that areas of radius 6° centred on the ten most over-dense pixels in Fig. are 
removed (sec table|2} corresponding to a removal of 6.3 per cent of the galaxies and 2.6 per cent of the solid angle of the full sample. 
Third column: Here, we split the full sample into north (triangles) and south (stars). For clarity, we have displaced the southern angular 
scaling parameter datapoints to the right. 
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I {°) 


b (°) 


n (dcg ) 


Cluster 


1 


312.2 


30.0 


69.1 


Shapley supercluster 


2 


266.5 


-51.3 


54.4 


Horologium-Reticulum 


3 


45.0 


57.4 


52.4 


Centaurus 


4 


10.4 


51.3 


52.2 


Centaurus 


5 


303.8 


32.8 


49.4 


Shapley supercluster 


6 


266.8 


-48.1 


47.5 


Horologium-Reticulum 


7 


219.4 


-32.8 


47.5 


NGC 1600 Group 


8 


343.1 


-32.8 


46.8 


Pavo-Indus wall 


9 


312.2 


35.7 


46.2 


Shapley supercluster 


10 


9.6 


48.1 


46.1 


Centaurus 



Table 3. The positions and galaxy densities of the ten most over- 
dense pixels of the smoothed Ks < 13.5 2MASS galaxy distribu- 
tion shown in Fig. |S] We also note the cluster with which each 
pixel is associated. 



|20o2); approximately, the scales used represent fits in the 
non-linear and quasi-linear or linear regimes. Also, it is im- 
portant to remember that the conversion from angular to 
real space scaling parameters becomes increasingly uncer- 
tain at large angular scales {6 >2°). The best fit scaling pa- 
rameters are shown in Table for small scales and Table |21 
for large scales. 

4.2 Fair sample issues 

It was noted in section 1, that a considerable problem 
in previous high-order clustering analyses is the presence 
of extreme fluctuations in the galaxy samples which have 
a significant effect on the observed scaling parameters. 
The APM and 2dF Galaxy Redshift Surveys observe ris- 
ing Sp values for r>4/i~^Mpc for example, consistent 
with some models of structu re formation with st rongly 
non-Gaussian inital consitions jGaztafiaga fc Fosalbp lTooi 
iGaztaiiaga fc Mahonenlll996l : lBemEgdeai^t aL.200^ : with 
the 2dFGRS at least the clustering signal is significantly al- 
tered when two superclusters are removed from the sample 
(corresponding to a reduction in the volume of ~2 per cent) 
yielding consta nt Sp parameters for p ^ 5 to scal es of 
r «10/i"^Mpc jBaueh et alJl2004l: ICroton et alJl2004^ . 

We therefore wish to examine the level of such sys- 
tematic effects by removing large over-densities from the 
2MASS galaxy sample and examining the effect on the ob- 
served scaling parameters. In Fig. |3 we show the pixelated 
2MASS Ka <13.5 galaxy density distribution smoothed on 
Rs7° scales (each pixel has a solid angle of 13.5 deg^); clus- 
ters are indicated by lighter filled contours and regions of 
under-density by darker shades. In the upper panel we show 
the entire Ks <13.5 2MASS sample in Aitofi^ projection; in 
the lower panels the galaxy density field for each galactic 
hemisphere is plotted in projection separately. We also list 
the ten most over-dense pixels (corresponding to a limit in 
the number density of n >46 deg~^) in table |21 

In Fig. |5| we plot [1)2 and the high-order angular scaling 
parameters with various cuts to the full galaxy sample used 
in Fig. In each case, the corresponding result for the full 
sample is indicated by a solid line. First, we omit the largest 
supercluster only; in the first column, we have removed a 
circular region with an angular radius of 6° centred on the 
most over-dense pixel in Fig. |3 sampling the Shapley su- 



percluster. This corresponds to a removal of 1.1 per cent of 
the galaxies and 0.33 per cent of the solid angle of the full 
|6| ^10° sample. The effect of this removal on the 2-point 
correlation function is insignificant. The form of S3 remains 
consistent with the result for the full sample, although the 
best fit slope at large scales (1°.0 < 8 <10°) changes to 
7 = -0.27toll (compared to 7 = 0.0llE!:^3 for the fuU 
sample). At higher orders, there is also no significant effect 
although the statistical uncertainty increases at large scales 
(datapoints with extremely large errorbars are omitted for 
clarity) . 

In the second column in Fig. |21 we omit all galaxies 
within 6° of the ten most over-dense pixel centres (see ta- 
ble l^l. This corresponds to a removal of 6.3 per cent of the 
galaxies and 2.6 per cent of the solid angle of the full sample. 
There is a small effect on the 2-point correlation function at 
large scales {9> 10°). The effect on the higher order angular 
scaling parameters also becomes more significant. However 
the effect on S3 at least is limited; the best fit slope at large 
scales (1°.0 < 9 <10°) is 7 = -0.33+o;2g. 

It has previously been observed that the form of 
the galaxy distribution in the northern and southern 
galactic caps are s ignificantly different in many respects; 
iMaller et all ll2005l) detected a difference in the b ^20° 
and 6 ^-20° 2MASS 2-point angular correlation func- 
tions at large scales for example. It has also recently been 
suggested that the southern galactic cap may contain a 
large local under-density in the galaxy distribution cov- 
ering > 4000 deg^ around the Southern Galactic Pole to 
r ~300 Mpc, which may be at odds with th e form of the 
ACDM P(k) at large scales jFrith et alJl2005al . 12004'. '20031: 
Bussw cU ct al. 2004 ). It is therefore interesting to compare 
the form of high-order clustering statistics in the galactic 
caps. 

We plot the 2-point angular correlation functions and 
high-order scaling parameters for b J5l0° and b ^-10° 
Ks <13.5 2MASS galaxies in the third column of Fig. H 
We observe a similar dis crepancy in the 2-po int function as 
determined previously bv lMaller et alJ i2005ll with a steeper 
slope in the southern lj2 at large scales. However, the high- 
order scaling parameters in the northern and southern local 
galaxy distributions, for p ^4 at least, are consistent with 
each other and constant values over ^3 orders of magnitude 
of angular scales to S< 20° . At higher orders there are differ- 
ences between the two which increase with p. Whether these 
differences are due to real north-south differences in the lo- 
cal galaxy distribution, or simply due to other systematic 
effects is unclear. 



5 DISCUSSION 

In the previous section, high-order angular and real space hi- 
erarchical scaling parameters were determined to extremely 
large scales (r<40ft~^Mpc) from a sample of 650 745 
Ks <13.5 2MASS galaxies. Unlike previous high-order clus- 
tering analyses, the scaling relations for p< 4 are relatively 
robust even when the most prominent clusters are removed 
from the sample. 

As described in sections 1 and 2, high-order clustering 
statistics have significant implications for the nature of the 
primordial density distribution and also the way in which 
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Figure 3. Filled contours representing the Itr, 2a and 3o" confi- 
dence regions for the real space galaxy skewness slope and ampli- 
tude determined from fits (accounting for the covariance) to 
the 2MASS |6| ^10° Ks <13.5 parameter (as shown in Fig.EJ 
at large scales (4.0< f <40/i~^Mpc). The best fit parameters 
are A=4.0 and 7s=-0.02, where we model the skewness using 
= A rTS . The dashed line indicates the predicted constant 
form (i.e. 7s=0) for 53 in the case of a Gaussian distribution of 
primordial density fluctuations. 



strongly non-Gaussian Initial Conditions 




0.0 0.5 1.0 1.5 

A 

Figure 4. Filled contours representing the Icr, 2a and 3(t confi- 
dence regions for the real space parameter (see eauation l2UI 
slope and amplitude determined from fits as in Fig. |3] The 
best fit parameters are yl=0.53 and 73 =-0.93, where we use 
B3 = A rTs as before. The dashed line indicates the predicted 
constant form for in the case of a strongly non-Gaussian dis- 
tribution of primordial density fluctuations. 



galaxies trace the underlying mass distribution. In the fol- 
lowing section we examine the consistency of our results with 
Gaussian initial conditions and hierarchical clustering and 
discuss the implications for non-Gaussian models of the pri- 
mordial density distribution. Assuming the former, we also 
infer constraints on non-linear galaxy bias. 

5.1 Implications for tlie primordial density field 

5.1.1 Consistency with Gaussian initial conditions 

The nature of the primordial distribution of density fluc- 
tuations is predicted to be close to or exactl y Gaussian in 
stand a rd inflationar y mode ls et al. 1993; Gangui ct alJ 

| l994l: iLesgoureud \\m\ IWang fc Kamionkowski .20001 
lMaldacenall2002l: lAcauaviva et alJl2003^ . This occurs as a 
consequence of the slow-roll conditions on the inflation po- 
tential which require the potential energy of the fleld to dom- 
inate over the kinetic energy in order to produce a phase 
of accelerated expansion which lasts for a sufflciently long 
period of time. Deviations from Gaussianity under these as- 
sumptions depend on the inflationary model but are gen- 
erally extremely small; in recent power spectrum analyses 
for example it is assumed that the primordial density field 
has ex actly random phases (e.g. lCole et alj2005l : lFrith et alJ 

Under the assumption of Gaussian initial conditions, 



a hierarchical scaling (see equation IIUI of the high-order 
moments of the dark matter and galaxy density fields 
iFrv fc Gaztanagal Il993h is expected through the evolu- 
tion under gravitational inst ability of the p rimordial den- 
sity fluctuations (e.g. Peebles 1980; Fry 1984: Bouche t et al] 
. 1992; Bcrnardcau 1992, 1994a; Gaztafiaga fc Baughl ll99a 
baugh et alJl2004h . This holds on scales where the gravita- 
tional collapse of dark matter haloes evolves linearly; it is 
important therefore to examine high-order moments in the 
linear and quasi-linear regime. 

It is clear from Fig. that our results provide an 
excellent agreement with the expected scaling relations 
at large scales, with Sp and Sp parameters consistent 
with constant values over several orders of magnitude 
in scale for p ^7. As noted previously by the 2dFGRS 
iBaugh et alJ l2004^ . this hierarchical scaling extends to 
small er scales than expected b y perturbation theory predic- 
tions (iBernardeau et alJl2002ft . This consistency with Gaus- 
sian initial conditions is also apparent from Tables and |5| 
where we perform fits to the scaling parameters to unprece- 
dented scales; there is good agreement between the small 
and large scale fits (although there is a small correlation be- 
tween the two). Alternatively, we can examine the measured 
slope of the Sp parameter at the scales of interest; in Fig.|3 
we determine confidence limits for the slope and amplitude 
of the S3 parameter at large scales (4.0< f <40/i~^Mpc) 
as shown in Fig. The constraints provide excellent agree- 
ment with the expected constant value as a function of scale. 
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(where Sa oc r~' 



with a best fit slope of 7s 
marginalising over the normalisation). 

This consistency with hierarchical scaling represents 
a departure from recent constraints on hi gh-order correla- 
tion functions. The APM llGaztafiagalll994l) and Edinburgh- 
Durh am Southern Galaxy Catalogue iSzapudi fc Gaztanaeal 
Il998l) surveys observe significant upturns in the an- 
gular skewness at scales of 6> 1° (corr esponding to 
r>4 /i ~^ Mpc). Also, ana lysis of the 2dFGRS llCroton et alJ 
I2OO4I: iBaugh et all l2004l) indicates rising Sp parameters. 
The issue is complicated by the fact that these measure- 
ments are not independent as there is significant overlap 
between the observed survey fields. Coupled with this, the 
2dFGRS identify two superclusters which significantly alter 
their results; when removed there is good agreement with 
the expected hierarchical scaling. However, due to the sur- 
vey volume the scaling parameters are determined only be- 
low r «10/i~^ Mpc and do not probe the larger scales of 
interest here. 



5.1.2 Constraints on non-Gaussianity 

While the assumption of Gaussian initial conditions is ac- 
ceptable assuming the validity of the simplest infiationary 
models, there also exist alternative models of inflation in 
which deviations from a Gaussian form for the primordial 
density field are predicted. For example, it is possible to in- 
troduce non-Gaussianity if the scalar field driving inflation 
has more than one component (e.g. iKofman fc PogosvanI 
Ll988; Bcrnardcau & Uzan 2003); strongly non-Gaussian ini- 
tial conditions are also a feature of models in which the 
infiaton is not a slow ly rolling scalar field but a fast mov- 
ing ghost condensate iArkani-Hamed et al.l20o3) . Addition- 
ally, it is possible to construct models in which the initial 
conditions are non-Gaussian as a result of non-linear struc- 
tures, such as cosmic strings or global textures, within the 
primordial density field (e.g . iGaztanaga fc Mahonenl|l996t 
lAvelino et aIlll998HGaztafiaga fc Fosalballl998D" 

When considering the effect of non-Gaussian initial con- 
ditions on moments of the local galaxy density field, it is con- 
venient to consider the departure from primordial Gaussian- 
ity in two ways. Here, we consider the strongly non-Gaussian 
regime as might arise from cosmic strings or textures; weak 
departures from Gaussianity as might arise from the vari- 
ous models of infiation described above are more difficult 
to constrain since the effect can be a shift in the Sp and 
Sp parameter amplitudes rather than a change in slope. In 
the case of strong departures from Gaussian initial condi- 
tions, the scaling of high-order moments of the density field 
is expected: 

^P = Bp (20) 

where for non- Gaussian models seeded by topological de- 
fects the Bp scaling parameters (not to be confused 
with Bp used in equation I15|l are exp ected to be con- 
stant at large scales and of order unity (jTurol^^S^e^ell 



ri99i]:lGazta,fiaga fc Mahonenll99(1:lRerna,rdeR,u et alJ200 

Equivalently, the typical signature of strong non-Gaussianity 
in the Sp scaling parameters is a slope of 7s>0.6 (where 
at large scales, with a characteristic minimum at 



B3 parameter through fits to the 2MASS results (mea- 
suring the projected -B3 = uj^./u^^'^ and using the transfor- 
mation to real space described in eguation llSH at large scales 
(4.0< f <40/i~^ Mpc). We are able to reject the scaling ex- 
pected in strongly non-Gaussian models (i.e. constant Bp 
parameters, 7s =0) as described in equation 1201 and there- 
fore non- Gaussian models seeded by topological defects, at 
the ~2.5a confidence level. 



5.2 Non-linear galaxy bias 

High-order clustering analysis also represents a powerful 
probe of the way in which galaxies trace the underlying mass 
distribution. We have previously noted that the galaxy bias 
associated with the variance of the density field has been 
constrained to h\ > 1 in the isTs-band at the > 3(7 level 
ijFrith et al.ll2005bl) and foi « 1 for optically-selected galax- 
ies. It is unclear whether deviations from the linear bias 
model exist, and if they do, how this coincides with current 
theories of galaxy evolution. 

Assuming Gaussian initial conditions, we use predic- 
tions from perturbation theory for the dark matter skewness 
and the relationship between the dark matter and galaxy 
skewness (equations 1121 and 113^ to compute the A"s-band 
non- linear bias coefficient C2 ~ &2/bi i mplied by our results . 
We use an n=-2 power spectrum slope ijPercival et alJlioOll : 
ICole et al.l2005h and a Ks-hwd linear bias of 6i=1.39±0.12 
measured from the Ks <13.5 2MASS galaxy angular power 
spectrum ijFrith et alJl200 5bV Since these predictions from 
perturbation t heory are valid only in the quasi-linear and 
linear regimes iBernardeaiJll994bl) . we use only the fits on 
large scales. First, we use the best fit galaxy skewness in 
the range 4.0< f <40/i"^Mpc (see Tabled; we find that 
C2=0.57±0.41. Since the conversion from S3 to S3 become s 
increasingly uncertain at scales oi 9 > 2° (|Gaztaiiagalll99 J ) . 
we also use a narrower range of 4.0< f <10h~^ Mpc (equiv- 
alent to 1°.0 < S < 2°. 5). Using the best fit galaxy skewness 
in this range of S'3=4.01±0.34, we find that C2 = 0.57±0.33. 

We therefore detect a positive quadratic contribution 
to the A"s-band galaxy bias at the ~ 2a level. This means 
that the density of 2MASS galaxies rises more quickly than 
the mass density contrast. This differs from all previous con- 
straints on the C2 non-linear bias parameter which have been 
negative, most recently with constraints from the 2dFGRS 
which limit bi = 0.94l° l? and C2 = -0.361^;^^ in the opti- 
cal 6j-band iGaztanaga et all2005h . and also from the IRAS 
PSCz catalogue o f bi = 0.83 ± 0.13 and C2 = -0.50 ± 0.48 
in the infra-red jFeldman et a l. 2001). This compares to 
near infra-r ed i^s- band constraints of 61 = 1.39 ± 0.12 
and the optimal constraint in this work 



Frith I 
of C2 = 0.57 ± 0.33. 

It is possible to understand these results by examining 



the analytic predictions of IMo et ahl H997.) for the high- 
order bias coefflcients, formed via the iPress fc SchechteJ 
il974fl formalism and an initially Gaussian density field. 
From this the first two terms in the Taylor expansion of 
eguation llll are predicted to be: 



bi = l + 



- 1 



ilOft^^Mpc ijGaztanaga fc Mahonenl ll996l. 
In Fig. 2] we constrain the slope and amplitude of the 



62 = 



(21) 
(22) 
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where f = 5c/(j{M), a{M) denotes the hnear rms fluctua- 
tion on the mass scale of a dark matter halo of mass M and 
Sc is the linear theory over-density at the time of collapse 
(for reference 5c=l-686 for Q—1). 

This relatively simplistic scenario is not able to pro- 
vide accurate quantitative predictions which match the ob- 
servational results above. However, by considering haloes of 
differing mass it is possible, qualitatively at least, to under- 
stand these apparently contradictory constraints on C2. For 
instan ce, if we conside r the most massive haloes for which 
iy^> 3. lMo et al.l lll997l) predict 61 >1 and C2 >0. In contrast, 
for typical mass haloes for which ~1, bias parameters of 
&i «1 and C2 <0 might be expected. With this in mind, it is 
possible to understand why a near infra-red survey, which is 
more sensitiv e to early-typ e galaxies than optical or infra-red 
surve ys (e.g. IJarrettlBool ICroton et alll2005l : lOliver et all 
Il996l . respectively), produces galaxy samples with higher 
values for the linear and quadratic bias parameters. 



6 CONCLUSIONS 

We have measured reduced angular correlation functions, 
ujp, to ninth order using 650 745 galaxies selected from the 
2MASS extended source catalogue. From our estimates for 
the angular correlation functions, we have determined the 
projected and real space hierarchical scaling parameters, Sp 
and Sp respectively. The prime motivation for such analysis 
is to test the hierarchical scaling hypothesis which predicts 
these parameters t o be constant in the linear and quasi- 
linea r regimes (e.g. |Peebleslll98d : lrHlll984l: iBouchet et all 
Il992l: lBernardeaulll992^ " As such we are able to probe the 
primordial density field and constrain various models of in- 
flation and structure formation. High-order clustering anal- 
ysis also allows us to probe the way in which galaxies trace 
the underlying mass distribution; a negative offset between 
ACDM predictions and observations by the 2dFGRS have 
recently been interpreted as evidence for a quadratic con- 
tribution to the galaxy bias, although these conclusions 
are based on con straints in the weakly non-linear regime 
jGaztaiiaga et al ]E005). 

The most com parable recent w ork are the analy ses of 
the 2dFGRS (Croton et al.^ 2004; Bau gh et akl |2004|) and 
APM surveys (Gaztafiaga 1994); the galaxy sample used 
in this work represents an order of magnitude increase in 
volume and solid angle over each respectively. Previous 
analyses of high-order clustering statistics have proved ex- 
tremely challenging; due to the relatively small volumes 
probed in 3-dimensional analyses for instance, the results 
have been sensitive to the presence of clusters and super- 
clusters within the galaxy sample. Not only this, but di- 
rect comparisons with perturbation theory have also proved 
difficult since the statistical uncertainty at large scales is 
considerable, and frequently the covariance in the statis- 
tics has been ignored leading to unrealistically small errors. 
In addition, the results for the hierarchical scaling param- 
eters have consistently displayed a puzzling upturn in val- 
ues at r ~4:h~^ Mpc, consistent with some models of struc- 
ture formation with strongly n o n-Gaussian initial conditions 
jGaztafiaga fc Mahonen||l996|: [ Gaztafiaga fc Fosalballl998l : 
IWhitelll999l:lBernardeau et al.ll2002h ! 

Here, we are able to determine the scaling parameters to 



high accuracy to unprecedented scales, r< 100 Mpc. We 
also carry out a full covariance analysis in order to take ac- 
count of correlations in the datapoints at different cell radii. 
We are therefore in a position to make direct comparisons 
with the predictions of perturbation theory since we probe 
well into the quasi-linear and linear regimes and we have a 
good understanding of the statistical uncertainty. However, 
since we are working with a projected galaxy sample, we 
have to convert the associated angular scaling parameters to 
real space via a transformation which becomes increasingly 
uncertain on large scales {9 >2°). Nevertheless, this work 
currently provides the best estimates of high-order cluster- 
ing statistics at large scales. We are able to reach a number 
of conclusions: 

(i) Our results are in line with the expected hierar- 
chical scaling relation, with Sp and Sp parameters con- 
sistent with constant values over ~3 orders of magnitude 
in scale to r ~40/i~^Mpc for p ^7; we constrain the 
slope of 53 to 7s = -0.02lgj* (where 5"^ oc r^). Such 
a scaling pattern is expected if an initially Gaussian den- 
sity field evolves under the action of gravitational insta- 
bilitv fe.g. jPeeble s 1980; Frv 1984: Bouchc t et aL .1991 
lBernardeaulll992lf . This result is in contrast to recent re- 
sults drawn frorn the 2dFGRS, APM and EDSGC surveys 
iGaztafiagalll994l: ISzapudi fc Gaztafiag3ll998l : lBaugh etal] 
I2OO4I : ICroton et alJ \2QQ4) which displav rising scaling pa- 
rameters at large scales. 

(ii) The scaling parameters are relatively robust to the 
removal of the largest over-density, the Shapley superclus- 
ter, although the best fit slope of S3 at large scales (1°.0 < 
6 <10°) becomes steeper, yet is still consistent with a con- 
stant value (7 = — O.27I0 29). We also use a more stringent 
cut by removing the ten most over-dense regions (see Fig. |S] 
and Table 13 corresponding to a cut to the main sample of 
6.3 per cent of the galaxies and 2.6 per cent of the total solid 
angle; the resulting S3 parameter remains broadly consistent 
with the result for the main sample, although with a slightly 
steeper best fit slope at large scales (7 — — 0.33lo 29). 

(iii) Since strong non-Gaussianity in the primordial 
density field, as might be expected in models seeded 
by topol ogical defects such as cosmic strings or global 
textures (|Avelino_et_al ] |l998l : iGaztafiaga fc Fosalbal Il99i : 
iGaztafjag^^Mahpnenll^g^ . results in a strong upturn 
in the hierarch ical scaling parameters at large scales (e.g. 
iBernardeau et al. 2002'). we are able to reject strongly non- 
Gaussian initial conditions, producing relations of the form 

^ ^2''^! 8,t the ~2.5o" confidence level. 

(iv) We compare our constraints on 5*3 at large scales 
(where we have assumed a constant value) to predictions 
from perturbation theory. We detect a significant deviation 
consistent with a non-linear, quadratic contribution to the 
Ks-hand galaxy bias, parameterised as C2=0.57±0.33 (de- 
rived from fits in the range 4.0< f <10/i~^Mpc), imply- 
ing that the 2MASS galaxy density rises more quickly than 
the mass density contrast. This positive result represents a 
significant difference from the negative values found previ- 
ously; constraints on C2 from the optically-selected 2dFGRS 
and the infra-red PSCz samples yield negative values. We 
explain these apparently con tradictory results through an 
examination of the model of iMo et al.l (^93), which pre- 
dicts bias parameters of &i > 1 (the linear bias) and C2 >0 
(the quadratic bias) if the surveyed galaxies typically reside 
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in large mass haloes. Similarly, as observed previously, we 
might expect bi ~1 and C2 <0 from galaxy samples in which 
late-type galaxies are over-represented. 
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Figure 5. The number density (in deg~^) of 2MASS Ks <13.5 galaxies binned in 13.5 deg^ pixels; under-dense regions are indicated 
by the dark filled contours; areas of over-density by lighter filled contours. For reference, the mean |6| ^ 10° number density is 19.1 
deg~^. In the upper plot we show the entire sky in Altos' projection; in the lower panels we show the b ^0° (left-hand plot) and b ^0° 
(right-hand plot) hemispheres, such that the galactic poles are positioned in the centres of each. Prominent clusters are indicated as 
follows: Co - Coma cluster, S - Shapley supercluster, V - Virgo supercluster, UM - Ursa Major cloud, GW - Great Wall, Ce - Centaurus, 
HR - Horologium- Reticulum, PP - Perseus-Pisces chain, PI - Pavo-Indus wall, N - NGC 1600 Group, Pu - Puppis, EF - Eridanus-Fornax. 



